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Abstract 

The graph of neutrinoless double beta decay is applied to HERA and generalized 
to final states with any two charged leptons. Considered is the case in which one 
of the two escapes typical identification criteria and the case when a produced tau 
decays hadronically. Both possibilities give one isolated lepton with high transverse 
momentum, hadronic activity and an imbalance in transverse momentum. We ex- 
amine the kinematical properties of these events and compare them with the high 
Pt isolated leptons reported by the HI collaboration. Their positive charged muon 
events can be explained by the "double beta" process and we discuss possibilities for 
the precise determination which original final state produced the single isolated lep- 
ton. To confirm our hypothesis one should search in the data for high pseudorapidity 
and/or low pr leptons or for additional separated jets. 
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1 Introduction 



Despite the impressive confirmation of predictions from the standard model (SM) it is 
general believe that we are on the verge of fundamental new discoveries, be it production 
of new particles or significant deviations of observables in high-precision measurements. 
Effects of new physics might also be hidden in existing data sets and it is interesting to see 
what candidates are able to explain any unexpected events or measurements. A first step 
in this direction was done in terms of observation of nonvanishing neutrino rest masses, 
most clearly seen in the up-down asymmetry of the atmospheric muon neutrino flux in Su- 
per Kamiokande UJ. The smallness of these masses can be related to massive new particles 
via the see-saw-mechanism 0. In this respect it seems most natural to look for effects 
of massive neutrinos, i.e. search for hints of these new particles in high- or low-energy 
experiments. The theoretical prejudice is that the neutrinos are Majorana particles — be 
it because they are delivered by see-saw or pop out of almost every GUT — and we shall 
follow this idea. 

Especially for the case of heavy (few 100 GeV) Majorana neutrinos, production at acceler- 
ators has been investigated by many authors. The different possibilities include e + e~ [[|, 
PP @; PP 0) U N 0, ep 0) HI) linear colliders and e~e~ |], |T(J or even e/i machines [[11[ . 



Heavy Majoranas have also been studied within the context of low-energy experiments 



such as neutrinoless double beta decay (0vf3(3) |)], [12] or Kaon decays JT3[. The respective 
Feynman diagrams as well as the concrete model differ in most publications and the inter- 
ested reader might compare the papers with respect to that. 

One of the anomalies in existing data is the existence of high px isolated leptons together 
with large missing transverse momentum (pj) at HERA. Since the first event [14] was 
discovered by HI, five more were found (T^j and at least 3 of them can not be explained 
by W production or other SM processes. In contrast to that, ZEUS sees no excess in these 



events [jT6| , yet, at the present statistical level, there is no contradiction [[17 
In H we examined the process (see Fig. 

e + p — > Til a + (3 + X with a, (3 = e, fi, r (1) 

and discussed possible signals of this like-sign dileptons (LSD) and high final state. 
No such events are reported and previously unavailable direct limits on the elements of 
the Majorana mass matrix were derived pi [I8| . However, it turns out that when the 
kinematical cuts used in Hi's search for isolated leptons are applied to our process (|l|), 
they tend to ignore one of the two leptons. Especially the requirement of p^ pton > 10 GeV 
is often too much for both charged leptons to fulfill. The LSD signal of Eq. ([!]) is thus 
reduced to one isolated lepton with high p T . This possibility can be checked by looking 
for an additional isolated low p? and/or high pseudorapidity lepton. In addition, it is 
possible that a produced r decays hadronically[], resulting also in single lepton final states. 
More than one isolated jet would be a signal for this kind of event. Since process ([I]) gives 
LSD with the same sign as the incoming lepton we concentrate on Hi's positive charged 

■""We shall use the term electron, muon or tau for both, particle and antiparticle. 
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muon events, since there are no positron events found. This fact might be explained if one 
incorporates also limits on mixing of heavy neutrinos, as derived from 0v/3(3. With this 
constraint the expected e signal is smaller than the /i signal. 

One might argue that direct Majorana (N) production via a e + NW vertex is more likely 
to occur since the cross section is larger. At present there is only an analysis in HERA's 
e~p mode available and it was found that detection is only possible if the N decays 
into e + W~, giving an isolated lepton with different charge than the incoming one. The 
reason for that is of course the large background from W production. However, a general 
analysis of all channels (N — > uZ, N — > fiW, . . .) remains still to be done and it might 
be interesting to compare the results with our signals in the future. Until that is done, 
we think that our process is worth considering, inasmuch as we have no restriction to the 
flavor of the final state lepton. 

The paper is organized as follows: In Section || we discuss some general features of the 
process and the diagram and argue in Section |3| that the two lepton signal of Eq. ([[]) might 
very well be seen as a one lepton signal. Section f| sees a discussion of signals of the events 
and how one might distinguish the genuine final state from the measured one. Finally, 
Section [5] closes the paper with a conclusion and discussion. 



2 The process and heavy Majorana neutrinos 

We shall work in a mild extension of the SM with no further specification of how heavy 
Majoranas might be created. The coupling to the usual leptons and gauge bosons is 
the familiar left-handed weak interaction. The three known light neutrinos u a are thus 
mixtures of light and heavy mass particles, this can be expressed by the replacement 

u a — > cos 8 a u a + sin 6 a N a (2) 

in the (unmixed) Lagrangian for each family. For the sake of simplicity we take N a = N. 
The Lagrangian now reads: 

- C = -£=W, t {coster 7 _Z Q + sin0 o JVy*7_Z a } 

v 2 (3) 

+ 2co'sfl Z » { cos2 e cP^fl- v a + sin 26 a NYl-Va - \ sin 2 9 a Nj^ 5 N} + h.c. 

where 7_ = |(1 — 75) and there is no vector current between iV and iV due to their 
Majorana nature. We can keep it for the light neutrinos since for energies much larger 
than the (light) masses there is hardly a chance to find a difference between 2V r )^ r y_v and 

777^75*/ m. 

What can we expect for the values of the masses and the mixing parameters? Taking the 
typical see-saw formula we find 



m 2 D (10 5 ...10 n ^ 2 

— - m N ~ f 

m N 10" 5 ... 1 



m u ~ —± m N ~ ^ ^ '- eV ~ 100 . . . 10 iS GeV (4) 
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where we took for the Dirac mass mj) every value from electron to top mass and for the 
light neutrino mass we allowed everything from the vacuum solution in a highly hierarchical 
scheme (a/ Am 2 ~ m v ~ 1CT 5 eV) to a degenerate scheme (cosmological or also LSND's 
mass scale) J2 m u — few eV (see |1| for a detailed analysis of allowed schemes). For the 
mixing angle we have 



6L ~ sin o 



m D 
M N 



10" 5 . . . 10 



-16 



(5) 



However, we shall use the current bounds on 9 a which are [20 



sin 2 9 e < 6.6 • 10~ J , sin 2 6„ < 6.0 • 10" d and sin 2 9 T < 1.8 • 10 



1-2 



(6) 



Note that the lowest value is for the muon sector. Eq. ([|) can now be applied to calculate 
the width of the Majorana, which is dominated by the two-body decays N — > Wa and 
N — > Zu n , we find 



r(A0 = £ 



Gp sin 2 C 
8nV2M 3 



N 



{[(il4 - M&) + M*,(M 



M 2 W ) 



(m; 



M 2 W ) 



+ cos o 



{M% - Ml) + Ml{Ml - Ml)] {Ml -Ml)}. 



(7) 



Direct searches for heavy neutrinos give typical lower limits |21| on their mass of 70 to 100 
GeV, depending on their character (Dirac in general gives a higher bound) and to which 
lepton family they couple to. Unfortunately, the maximal value of the cross section of 
process ([]]) in Fig. [1] is found to lie in that range as well 0, §. The dependence on the 
mass goes as 

for M 2 N < q 2 
l n 2 for M 2 N > q 2 

where q is the momentum of the Majorana. The standard calculation gives for the matrix 
element || (see Fig. |l| for the attachment of momenta): 



da oc 



M 2 N 
(q 2 - M] 



If 



Ml 



(8) 



\M\ 2 \ 



e^q 



v e a + a + q') 



s\v A 9 a M 2 N G 4 F Ml2 12 j^ 



1 

7TTT7ZT 



MwY(<& 



l WJ 



(4 ~ M\ 



I? 



(k 2 -Pi){h ■ h) + 



(a 



1 



M 



TT2~ 

NJ 



(h -Pi)(k 2 ■ k 4 



(9) 



(q 2 -M 2 N )(q~ 2 



M 



tt((^2 • h)(pi ■ h) - (h -pi)(h ■ h) - (k 3 -p 1 )(k 2 ■ fc 4 )) 
n) 



and the scattering with an antiquark sees k^ interchanged with p 2 . Here q 2 denotes the 
momentum of the Majorana in the crossed diagram, which has a relative sign due to the 
interchange of two identical fermion lines. In addition one has to include a factor | to 
avoid double counting in the phase space integration. For the phase space we called the 
routine GENBOD |2l| and for the parton distributions we applied GRV 98 [2~3||. In case a 
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r is produced we additionally folded in its three-body decay. We inserted finite (W and 
N) width effects in our program and found them to be negligible. 

An interesting statistical effect occurs when one considers the relative difference between, 
say, the and the [it final state (mass effects play no significant role): First, there is 
no factor \ for the latter case. Then, there is the possibility that a r is produced at the 
( "upper" ) e + T7^ W vertex or at the ( "lower" ) qq'W vertex. Both diagrams are topologically 
distinct and thus have to be treated separately. This means, four diagrams lead to the 
[it final state, whereas only two lead to the [i[i final state. We see that there is a relative 
factor 4 between the two cases. Note though that now the interference terms are added to 
the two squared amplitudes since there is no relative sign between the two. This reduces 
the relative factor to about 3. However, effects of kinematical cuts and the limits of Eq. 
(0) wash out this phenomenon. A similar situation occurs when one studies the rr case 
and lets the r's decay into different particles (e.g. evv and [wv). There is no way to tell 
into what the "upper" or "lower" tau decays, so one has to include both cases. 
A question arises if one can conclude a Majorana mass term if we measure a process like 
Eq. (|IJ). Here, a simple generalization of arguments first given by Schechter and Valle 
for neutrinoless double beta decay applies: Assuming we found indubitable evidence for 
e + q — > zTg a + /3 + q', then crossing permits the process — > e~Vl a + (3 + q'q, realized by the 
"black box" in Fig. |^. Any reasonable gauge theory will have W's couple to quarks and 
leptons, so that a Majorana mass term for u a and i/g is produced by coupling one W to 
the a + e~z7e and one W to the (5 + q'q vertex. Since we do not know which two quarks 
participate and which neutrino couples to the positron (the Schechter and Valle argument 
for works with two pairs of u and d quarks), this theorem holds for a greater class 

of models, namely e.g. those with direct e~ Vx coupling, with X being any flavor. 
The connection between a neutrinoless double beta decay signal and Majorana masses has 
been expanded in |25] to supersymmetric theories and it was found that it implies Majo- 
rana masses also for sneutrinos, the scalar superpartners of the neutrinos. 
However, in contrast to the signal in neutrinoless double beta decay experiments (two 
electrons with constant sum in energy) the identification of our process will be very dif- 
ficult and deciding which a(3 final state was originally produced remains a hard task. In 
addition, the number of expected events turns out to be far less than one. Nevertheless, 
the demonstration of Majorana mass terms will be an exciting and important result, since 
different models predict different texture zeros in the mass matrix. In some models the ee 
entry in the mass matrix is zero and therefore the only direct information about the mass 
matrix might come from neutrino oscillations, cosmological considerations, global fits and 
direct searches, e.g. at LEP. This complicates the situation, since e.g. in oscillations only 
mass squared differences are measured and the additional phases induced by the Majorana 
nature are unobservable. To combine all information from the different approaches input 
from models is required. Thus, the exact form of the matrix is highly nontrivial to find. 
Therefore, information about non-vanishing entries in the mass matrix is very important 
and one has to take every opportunity to find out about all elements and the Majorana 
character in general. In addition, if such a lepton-number violating process is detected, it 
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is surely helpful to know if the "mildly extended" SM can provide the signal or another 
theory, such as SUSY, has to be blamed. 

3 Two become one 

We applied the same cuts as HI in their search for isolated leptons: 

• Imbalance in transverse momentum jpj, > 25 GeV 

• Transverse momentum of lepton px > 10 GeV 

• Pseudorapidity of lepton < 2.436 

• Distance between charged lepton and closest jet in rj-cf) spaceQ, Ai? > 1.5 where (f) is 
the azimuthal angle 

• Angle of the hadronic jet(s) 4° < 6 X < 178° 

This has to be compared with our cuts in ||, pj, > 10 GeV, \rj\ < 2.0 for all measured 
particles and AR > 0.5 between the charged leptons and the hadronic remnants. It turns 
out that both sets of cuts deliver cross sections in the same ballpark. There are now two 
possibilities for the original LSD signal to appear as one single lepton: 

1. One lepton can have high pseudorapidity and/or low pt- This lepton then also 
contributes to the missing transverse momentum. 

2. If one tau is produced it might decay hadronically, adding one neutrino to the im- 
balance in px and also additional hadronic jets. 

A detailed analysis of the LSD signal might be done if one finds such events. Collecting 
all possibilities for the otf3 final states and the r decays results in Fig. ||. We denote with 
"hadronic" the signal coming from final states which also have additional hadronic activity 
from a r decay. We call "leptonic" the signals coming from events in which two final state 
leptons are produced from which one escapes the identification criteria. Those included 
therefore most channels, namely all of them except the ones with hadronic tau decay. One 
can see that muon events have a smaller cross section than the e signal, coming from the 
fact that their mixing with the heavy neutrino has the biggest constraint. Mass effects play 
no significant role. If the HI anomaly is indeed explained by heavy Majorana neutrinos, 
one might ask why only muon events are detected. A possible reason for that might lie 
in the following fact: The experimental constraint from Oz//5/5 on mixing with a heavy 
Majorana neutrino reads || 

sin 2 £ e < 5 ■ 10- 8 m_ (10) 

2 Actually Hi's value is 1.0 or 0.5, depending on the way they define jets for their respective analysis. 
We use a general value of 1.5 to account for hadronization effects. 
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Incorporating this in Fig. [3] gives Fig. [|. Now the electron signal is far below the muon 
signal. In Fig. [| we plot how the cross section for the production of a /x is composed. The 
biggest contribution comes from the /ir channel, which has its reason in the mentioned 
factor ~ 3 relative to the aa channels and the high hadronic tau branching ratio, BR(r — ► 
v T hadrons) ~ |. 

An additional reason why the ra channels add higher contributions is that when the tau 
decays it distributes its momentum to three particles, i.e. the pr is in general lower and it 
can therefore escape the pt > 10 GeV cut more easily and thus has a higher cross section. 
Addmittedly, the process gives only a tiny signal: Multiplying the cross section with the 
36.5 pb -1 luminosity HI analyzed, gives 10~ 8 . . . 10~ 9 . However, as explained at the end of 
the previous section, one has to check every possible appearance of Majorana mass terms 
in order to get information about the mass matrix. 

One sees that many LSD signals produce the same single lepton signature. In the next 
section we will discuss possibilities to distinguish different original final states from the 
measured ones. 



4 Signals and observables 

Some kinematic quantities of Hi's positive muon events are given in Table [l]. In their 
analysis ]T3| using 36.5 ± 1.1 pb -1 luminosity at E p = 820 GeV and E e = 27.5 GeV, 6 
events were found (0 e + , 1 e~, 2 2 and 1 fi of undetermined charge), where about 
2 e and 1 /i are expected from SM processes. From those, the most important ones are 
W production, NC events (for e + events) and photon-photon interactions (/i 1 * 1 ). We also 
include the event with undetermined charge. The e~ and one which also has a e + , 
are very likely to stem from W production. In the following we will plot all distributions 
for = 200 GeV since the qualitative conclusions we draw remain valid for all masses 
considered. In our analysis it turned out that — with the kinematics from Table |I| — scatter 
plots of pt, jp? and the transverse mass Mt = J (j>p + Pt) 2 — {fir + Pt) 2 are most useful. 
We stress again that for different final states the composition of the missing transverse 
momentum can be made of two particles (for fifi or ee final states) to 6 (1 lepton and 5 
neutrinos, rr channel with two leptonic decays) thus changing the area in which events 
populate, say, the Mt~Pt space. The transverse momentum is also very sensitive on the 
original final state since r decays share the initial momentum to three particles thereby 
reducing the average Pt- This is displayed in Fig. || (/i/i final state) and Fig. [7| (rr, one 
hadronic decay). If both taus decay leptonically the distribution looks similar. Obviously, 
the rr case has in general low px and Mt, whereas the \x\x case displays an uniform 
distribution with a slight band in the center region indicated. 

Turning now to pj, we see in Figs. [| (/^*) and § (/ir, hadronic decay) that the situation is 
not as clear in the "mixed" channels, i.e. channels with two different final state charged 
leptons. Though the population in p t -Mt space is different (lower values in the latter 
case), it is not as obvious as for the /i/i/rr case. Due to its low M T , event /ii seems to be 
less probable in all figures, though no definite statement can be made. 
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Hi 


H2 


A*5 


Pt 


23.4i£j? 


28.0i«;I 


> 44 




18.9t^ 


43.21^ 


> 18 


M T 


3.o±ss 


99 Q+^.Y 
ZZ.O_4 2 


> 54 


8 


18.9^ 


i7.ii?:? 


> 22 


Pt 


42.2 ±3.8 


67.4 ±5.4 


30.0 ±3.0 



Table 1: Kinematical quantities of Hi's candidates for positively charged muon events. 
Note that the charge for /i 5 is undetermined. All values in GeV, taken from \TE\. The 
limits for /i 5 correspond to 95 % C.L. 



Now we consider the quantities connected with the hadronic remnants. We found that 
the muon signal is composed of roughly 1/3 purely leptonic final states and 2/3 events with 
additional jets. An interesting quantity is 5 = E%{1 — cos 6i) where the sum goes over all 
measured final state particles. In Figs. |lfj and [□] one sees that the presence of three jets 
keeps 5 in more or less in the same area whereas the transverse momentum of the hadronic 
system p?f is shifted towards lower values. Here also fii hes in a less crowded area. 
What are now the signals of the escaping charged lepton (if there is one)? In Figs. O and 



13| we plot the pseudorapidity r\ of the undetected lepton against its transverse momentum. 



Again, the /z/z and the rr case can be distinguished since the latter has far lower p?. The r 
boosts its decay products more or less in its forward direction so that r\ does not alter much. 
The problem one might encounter is that the escaping lepton is hiding in the hadronic jet. 
Demanding a distance of AR > 1.5 reduces the cross section by 10 to 15 % but does not 
change the distributions in Figs. [T2| and |T3[ 

We mimicked the hadronic r decay via two quark jets and ignored effects of modes like 
t — > 7r's v t . Thus, due to the boost of the r, two of the three jets of events with a 
hadronically decaying tau will be very close together. Fig. T\ displays the normalized 
distribution of the distance in r]-(p space. We denote with Ri the distance ordered with 
ascending value. One distance is centered significantly below one, therefore, probably two 
jets instead of three will be measured. However, the r identification is hard to do and 
Fig. O serves only as an indication of how things might work. Information on the jet 



multiplicity is not given in Ref. [15], though Hi an d (15 seem to have additional separated 



tracks in their event displays as can be seen in Fig. 2 of Ref. |T5 



5 Conclusions 

In the light of recent data we did a full analysis of the analogue of the neutrinoless double 
beta decay graph with all possible two charged leptons in the final state. One lepton can 
escape the identification criteria by either having low transverse momentum and/or high 
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pseudorapidity or (if it is a r) via hadronic decay. Signatures of these events are discussed 
and compared to Hi's isolated leptons with large missing transverse momentum. All their 
positive muon events lie in the regions typically populated by the "double beta" process, 
though fix is always in a less crowded area. Due to its high errors ^ is mostly in the favored 
region (but for the same reason always in the region populated by SM processes ||15||). 
Though the process represents an attractive explanation, the smallness of the expected 
signal might spoil our interpretation. Nevertheless, any information about mass matrix 
entries and Majorana particles is very important and worth looking for, regardless of the 
small expected signal. Other extensions of the SM might give larger signals to the discussed 
final states and it is then helpful to know how the "SM + heavy Majorana" extension 
contributes. 

To confirm our hypothesis direct production of heavy Majorana neutrinos will be the only 
possibility since cross sections or decay widths of other 0^/3/5 like processes are probably 
too small to be detected RTH1- Here, either HERA itself or LHC will be the candidates for 



this observation. We did not consider the mass reconstruction of the Majorana since the 
large number of unmeasured particles does not permit that. 

Other proposed explanations Jl4| for the events were FCNC interactions (topologically 
identical to leptoquark production) or high jets from which one fakes a muon signal. 
Production of supersymmetric particles is suggested in p6[ to explain the results. A definite 
answer regarding all detector/identification issues can only be given by the collaboration 
itself. From the "new physics" side we believe that massive neutrinos provide one of the 
most natural possibilities. 
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Figure 1: Diagram for e + p — > v e a + f3 + X. Note that there is a crossed term and for a ^ j3 
there are two possibilities for the leptons to be emitted from. 
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Figure 3: Cross section of the expected isolated lepton signal from all possible final states 
as a function of the Majorana mass mjv. 
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Figure 4: Same as previous figure if one also incorporates the Oz//5/3 limit from Eq. (FLO 
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Figure 5: Contribution of the different final states to the cross section of the muon signal 
as a function of the Majorana mass m^. 
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Figure 6: Scatter plot of p T of the measured lepton against M T for the channel e + p 
V~ e n + ji + X with one muon escaping the identification criteria. 
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Figure 7: Scatter plot of px of the measured lepton against My for the channel e + p 
VI t + t + X with one tau decaying hadronically. 
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Figure 8: Scatter plot of p T against M T for the channel e + p — > z/ e fi + fi + X with one muon 
escaping the identification criteria. 
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Figure 9: Scatter plot of jp T against M T for the channel e + p — > v e r + fi + X with the tau 
decaying hadronically. 
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Figure 10: Scatter plot of p* against 5 for the channel e + p — > v e fi + fi + X with one muon 
escaping the identification criteria. 
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Figure 11: Scatter plot of against 5 for the channel e + p — > z/ e fi + r + X with the tau 
decaying hadronically. Note that the range is half as wide as in the previous figure. 
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Figure 12: Scatter plot of pt of the escaping charged lepton against its pseudorapidity for 
the channel e + p — > fi + fx + X. The sharp edge is due to the applied cuts from HI. 
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Figure 13: Scatter plot of pt of the escaping charged lepton against its pseudorapidity for 
the channel e + p — > VI t + t + X. The sharp edge is due to the applied cuts from HI. 
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Figure 14: Distance of the three jets in rj-cf) space with Ri ordered with ascending value 
for the process e + p — > V e r + T + X with one r decaying hadronically. 
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